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Abstract. In this paper we study the dominant microscopic processes that generate 
nearly the whole one-electron removal and addition spectral weight of the one- 
dimensional Hubbard model for all values of the on-site repulsion U. We find that 
for the doped Mott-Hubbard insulator there is a competition between the microscopic 
processes that generate the one-electron upper-Hubbard band spectral-weight 
distributions of the Mott-Hubbard insulating phase and finite-doping-concentration 
metallic phase, respectively. The spectral-weight distributions generated by the non- 
perturbative processes studied here are shown elsewhere to agree quantitatively for 
the whole momentum and energy bandwidth with the peak dispersions observed by 
angle-resolved photoelectron spectroscopy in quasi-one-dimensional compounds. 
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1. INTRODUCTION 

For energies larger than the transfer integrals for electronic hopping between the chains, 
the one-dimensional (ID) Hubbard Hamiltonian is the simplest model for the description 
of electronic correlation effects on the spectral properties of quasi- ID compounds 
PiaiSlll!. It reads, 

H = -^EH.-^i+i.- + ^-c-] + ^E%T%i , (1) 

where c] and Cj^a are spin-projection a =t, i electron operators at site j = 1, 2, A^^, 
%,o- ~ '^i <T Cj,o-, and t is the transfer integral. In contrast to other interacting models j3] 
and in spite of the model exact solution jEl [Z| , for finite values of the on-site repulsion U 
little is known about the non-perturbative microscopic processes that control its finite- 
energy spectral properties. Recently, the problem was studied in Refs. jHl E] by the use 
of a pseudofermion description. The preliminary predictions of the method introduced 
in these references agree quantitatively for the whole momentum and energy bandwidth 
with the peak dispersions observed by angle-resolved photoelectron spectroscopy in the 
quasi- ID conductor TTF-TCNQ [TJ |2] and are consistent with the phase diagram of 
the (TMTTF)2X and (TMTSF)2X series of compounds 4j. More recently, results for 
the TTF-TCNQ spectrum consistent with those of Refs. |2] were obtained by the 
dynamical density matrix renormalization group method [Sj. Within the method of Refs. 
|Hl E! , the finite-energy spectral properties are controlled by the overall pseudofermion 
phase shifts, through the pseudofermion anticommutators |ID|. Such a method is a 
generalization for all values of the on-site repulsion U of the technique introduced in 
Refs. muni for ?7 oo. 

The studies of this paper focus on the specific case of the one-electron spectral 
weight and use the pseudofermion description used in the general studies of Ref. [S]. 
Such a description is related to the holon, spinon, and c pseudoparticle representation 
introduced in Ref. |T3j. Here we study the generators of the one-electron dominant 
microscopic processes and the interplay between such processes and the metal - Mott- 
Hubbard insulator quantum phase transition ^Ij . Since the mechanisms found here are 
expected to occur in other correlated models, our results are of interest for the further 
general understanding of the microscopic mechanisms associated with the quantum 
phase transitions. Furthermore, they are of interest for the further understanding of 
the unusual spectral properties observed in low-dimensional materials. 

The model describes spin-up electrons and Ni spin-down electrons in a 
chain of Na sites. We denote the electronic number by N = + N^. The number 
of lattice sites Na is even and very large. For simplicity, we use units such that both 
the lattice spacing a and the Planck constant are one. In these units the chain length 
L is such that L = Naa = Na- Our results refer to periodic boundary conditions. 
We consider an electronic density n = + n-i in the range < n < 1 and a spin 
density m = — = 0, where ria = N^/L and a =t, |. We introduce the Fermi 
momenta which, in the thermodynamic limit L oo, are given by ^kpa = ±7mo- and 
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±kF = ^Ikp] + ^Fi]/2 = ±7m/2. The one-electron spectral function B\k, u) such that 
/ = — 1 (and / = +1) for electron removal (and addition) is given by, 

B-\k,uj)= \{f\ck,a\GS)\^6{L0 + Ef,r^^i-EGs), uj<0; 

B^\k,u;)= \{f'\clJGS)\' 5{lu - Ef,,r,+, + Egs) , uj > . (2) 

/' 

Here Ck^a (and c|. ^) are electron annihilation (and creation) operators of momentum k 
and \GS) denotes the initial A^-electron ground state. The / and /' summations run 
over the — 1 and + 1-electron excited states, respectively, and [i?/, Ar_i — Egs] and 
7V+1 — Egs] are the corresponding excitation energies. 

The Hamiltonian commutes with the generators of the ?7-spin and spin SU(2) 
algebras [131 US Hn]. Here we call the r^-spin and spin of the energy-eigenstates rj 
and S, respectively, and the corresponding projections rjz and Sz- All such states 
can be described in terms of occupancy configurations of r^-spin 1/2 holons, spin 1/2, 
spinous, and ?7-spin-less and spin-less cO pseudoparticles [13]. (The latter objects are 
called c pseudoparticles in Refs. [IHl IIZ|-) We use the notation ±1/2 holons and 
±1/2 spinous according to the values of the ?7-spin and spin projections, respectively. 
For large values oi U/t the ±1/2 holons and —1/2 holons become the holons and 
doublons, respectively, used in the studies of Ref. [TH]. The electron - rotated- 
electron unitary transformation pH! maps the electrons onto rotated electrons such 
that rotated-electron double occupation, unoccupancy, and spin-up and spin-down 
single occupation are good quantum numbers for all values of U. While the —1/2 
and ±1/2 holons refer to the rotated-electron doubly occupied and unoccupied sites, 
respectively, the —1/2 and ±1/2 spinous correspond to the spin degrees of freedom 
of the spin-down and spin-up rotated-electron singly occupied sites, respectively. The 
charge degrees of freedom of the latter sites are described by the spin-less and 77-spin-less 
cO pseudoparticles. In turn, the cz/ pseudoparticles (and sz/ pseudoparticles) such that 
V = 1,2,... are r^-spin singlet (and spin singlet) 2z/-holon (and 2z/-spinon) composite 
objects. Thus, the numbers of ±1/2 holons {a = c) and ±1/2 spinous {a = s) 
read Ma,±i/2 = -Z^a,±i/2 ± J2'^=i^Nau where a = c, s, Nau denotes the number of 
ap pseudoparticles, and I^c, ±1/2 = ± Vz] and i^s,±i/2 = ["S" ± S^] gives the number of 
±1/2 Yang holons and ±1/2 HL spinous, respectively. Those are the holons and spinous 
that are not part of composite pseudoparticles. The total number of holons (a = c) and 
spinous (a = s) is given by = + Mq,__i/2]. 

The cO pseudofermions and composite pseudofermions are generated from the cO 
pseudoparticles and composite pseudoparticles of Ref. [13] by a unitary transformation 
[Hj. It introduces shifts of order 1/L in the pseudoparticle discrete momentum values and 
leaves all other pseudoparticle properties invariant. It is useful for our study to consider 
the pseudofermion subspace (PS), which is spanned by the initial ground state \GS) 
and all excited energy eigenstates contained in the one-electron excitations Cj^^\GS) 
and Cj^„\GS) [Hj. The local au pseudofermion creation (and annihilation) operator 
fxj,au (and fxj,au) creates (and annihilates) a au pseudofermion at the effective au 
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lattice site of spatial coordinate Xj = j a^^u- Here j = 1,2, N*^^ and a°j, = l/n^^i, 
is the effective au lattice constant given in Eq. (55) of Ref. J7] in units of the 
electronic lattice constant and n*^ = N^^/L = N*^/Na. The general expression of 
the number of effective au lattice sites N*^^ is given in Eq. (B6) of Ref. |13jj, where 
the number of au pseudofermion holes N^^^ is provided in Eq. (B.ll) of the same 
reference. All PS energy eigenstates can be described by occupancy configurations of 
au pseudofermions, —1/2 Yang holons, and —1/2 HL spinous 13 . For the ground state, 
Nco = N, Nsi = Ni, and N^u = Ns^i = I^a,-i/2 = for a = c, s, z/ > 0, and u' > 0. 
The deviations AN^„ ANj^^, AL,,_i/2, AL,,_i/2, AM,,_i/2, AM,,_i/2, AM,, AM, 
of the above numbers which result from the ground-state - excited-energy-eigenstate 
transitions play a major role in our study. It follows from the values of the ground-state 
numbers that AA^ci. = Ncu, AN^u' = Nsu', ALc,_i/2 = -^^c,-i/2, ALs,„i/2 = -^^s,-i/2, 
AMc,_i/2 = Mc,_i/2, and AM,,_i/2 = M,,_i/2 for i/ > and z/' > 1. Thus, often 
we replace the latter deviations by the corresponding numbers. A concept widely used 
in our studies is that of a CPHS ensemble subspace jSUHl- (Here CPHS stands for cO 
pseudofermion, holon, and spinon.) Such a subspace is spanned by all energy eigenstates 
with fixed values for the —1/2 Yang holon number Lc, _i/2, —1/2 HL spinon number 
-^5,-1/2, and for the sets of au pseudofermion numbers {A*',,^} and {Ng^'} corresponding 
to the u = 0,1,2,... and u' = 1,2,... branches. Fortunately, nearly the whole one- 
electron weight corresponds to subspaces involving the cO, cl, si, and s2 branches only. 

For the Af = 1 electron problem, the operator 0{j- • = 0[ • of the spectral-function 
expressions of Eq. (7) of Ref. [Hj is the operator Cj o- for / = —1 and c]^ for / = +1. 
These expressions can be re-expressed in terms of the operator 6^ j, which plays an 
important role in our study and is defined in terms of the original AA-electron operator 
Oj^j- j in Eqs. (27) and (28) of Ref. [9 . For the present A/" = 1 problem we use the 
notations 9j^a and Oj^^ for the operators and Qfj, respectively, and call Oij^^j and 
Oj j^^ the operators 0]~|j and Qf^j, respectively, on the right-hand side of Eq. (32) of 
Ref. 0. The latter equation then reads, 

_ oo ^ 

=^0,i,a + EVC*^'^~i,i,-; J = 1,2,..., AT,; l = ±l, (3) 
1=1 

and a similar expression with c~^ replaced by cf^ holds for ^, where for i > the 
index i = 1,2,... is is a positive integer number which increases for increasing values 
of the number of extra pairs of creation and annihilation rotated-electron operators in 
the expressions of the operators Oi relative to that of Oqj^„ and the value of the 
constants cf^ reads Cq^ = 1 and for i > is a function of n, m, and U/t such that 
cf^ — > as ?7/t — *• oo p. Moreover, the operators 9j^„ and 6'o,j,cr of Eq. Q have the 
same expression in terms of local creation and annihilation electron and rotated-electron 
operators, respectively. 

According to the general studies of Refs. [1111], the i = operators do j^cr and Ooj^^ 
generate nearly the whole spectral weight of the corresponding one-electron problems 
for all values of U/t and L ^ oo. In the ensuing section we confirm numerically 
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that the same occurs for finite values of the chain length L and treat the problem 
analytically for L —>■ oo and large values of U/t. Since the physics associated with 
the dominant microscopic processes that control the one-electron weight distribution is 
simplest understood in terms of pseudofermion occupancies, in Sec. Ill we calculate the 
expressions of the above operators in terms of local pseudofermion operators. In Sec. IV 
we study the main microscopic effects of the metal - Mott-Hubbard insulator quantum 
phase transition onto the one-electron spectral properties. Finally, the concluding 
remarks are presented in Sec. V. 

2. ONE-ELECTRON SPECTRAL- WEIGHT DOMINANT PROCESSES 

Application onto the ground state of the operators Ol j^^ or 6'o,j,cr generates transitions 
to excited energy eigenstates whose Mc^_i/2 = i^c, -1/2 + J2'^=i^Ncu values obey the 
following exact charge selection rule [0], 

Mc, -1/2 = , electron removal ; _i/2 = 0,1, electron addition , (4) 

where electron refers here to rotated electron. Thus, in the present case we only need 
to consider excited states such that _i/2 < 1 and Ncu = for z/ > 1. Moreover, the 
values of the numbers of —1/2 HL spinous generated by application onto the ground 
state of the operators c] and cj^^ are restricted to the following ranges, 

Ls, -1/2 = , I electron removal, t electron addition 

= 0,1, I electron addition, | electron removal , (5) 

whereas the permitted values of i^c, -1/2 coincide with those of Eq. (jlj). 

While the selection rules (jH) and ^ are exact for the one-rotated-electron and 
one-electron problems, respectively, direct evaluation of the weights by the method 
introduced of Refs. jHl IH] reveals that for A'o = 1 94% to 98% of the spectral weight 
generated by the operators Oqj^^ or fi'o,^,^ corresponds to transitions to excited energy 
eigenstates whose i^s, -1/2 and Nsu values for i/ > 1 are in the following range, 

00 

Ls, -1/2 + X! ('^ ~ ^s'^ = ) i electron removal, | electron addition 

= 0,1, I electron addition, | electron removal , (6) 

where electron refers here to rotated electron. Thus, in the present case most of the 
spectral weight corresponds to excited states such that _i/2 < 1 and Ng^, = for 
u>2. 

One-electron processes associated with excited energy eigenstates whose deviations 
do not obey the rule (jH) are generated by the operators 9ij^a or Oj such that i > 0. 
We confirm below that application onto the ground state of all such i > operators 
amounts for less than 1% of the [k, u;)-plane one-electron removal or addition spectral 
weight. We start by confirming that the one-electron processes which generate excited 
energy eigenstates obeying the ground-state charge selection rule (j^ are dominant for 
finite values of L and correspond to 99.75% - 100.00% of the whole electronic spectral 
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weight for all values of U/t. For finite values of L the thermodynamic Bethe-ansatz 
equations introduced by Takahashi [7j do not apply, whereas the Bethe-ansatz solution 
and the concept of a rotated electron remain valid. Therefore, for finite values of L we 
can characterize the excited energy eigenstates in terms of rotated-electron occupancy 
configurations and thus of the related ±1/2 holon, ±1/2 spinon, and cO pseudofermion 
number deviations 17J. Moreover, for zero-magnetization initial ground states the 
number of excited-energy-eigenstate si pseudofermion holes equals the number of holes 
in the corresponding Bethe-ansatz spin excitation spectrum. Thus, one can classify the 
finite-L processes by the deviation values AA^^co, _i/2, AMs^_i/2, AMc, AMg, and 
AiVi\. 

The correlation energy 2fi, defined in Eq. (21) of Ref. jH], plays an important role 
in the model finite-energy spectrum. The limiting values of such a correlation energy are 
given in Eq. (13) of Ref. [0]. At half filling it equals the Mott-Hubbard gap, 2/i = Emh- 
Our considerations refer mostly to one-electron addition. For simplicity, let us assume 
that the initial ground state has zero spin density. In this case the spectral-weight 
distribution associated with creation of a spin-up electron has the same form as that 
associated with creation of a spin-down electron. Here we consider the former case. For 
creation of a spin-up rotated electron, the selection rule (jlj) only allows the two values 
AMc, _i/2 = Mc^_i/2 = 0, 1. From the relation of the electron and rotated-electron 
numbers to the holon, spinon, and cO pseudofermion numbers ^Sl ^] , we find that for 
creation of a spin-up rotated electron the following transitions obey the selection rule 
0: 

(i) Lower-Hubbard band (LHB) transitions such that AA^^co = 1, AMc^_i/2 = 0, 
AMs _i/2 = 0, AMc = —1, and AMg = 1. The minimal excitation energy for such 
transitions is zero. The sub-class of these transitions that also obey the restrictions of 
Eq. (0) are such that AA^^'\ = 1. 

(ii) Upper- Hubbard band (UHB) transitions such that ANco = —1, AM^ _i/2 = 1, 
AMs _i/2 = —1, AMc = 1) and AMg = —1. The minimal excitation energy for such 
transitions is 2/i. The sub-class of these transitions that also obey the restrictions of 
Eq. (P) correspond to a si pseudofermion-hole deviation value AA^^^ = 1. 

The simplest transitions that do not obey the selection rule involve creation of 
three cO pseudofermion holes and three holons: 

(iii) Second-UHB transitions such that AA^co = ^3, AMc^_i/2 = 2, AM^ _i/2 = —2, 
AMc = 3, and AMg = —3. The minimal excitation energy for such transitions is 4/i. 
The sub-class of these transitions that also obey the restrictions of Eq. © are such that 
AN^, = 1. 

The simplest transitions of types (i) and (ii) that do not obey the restrictions (jHI) 
involve creation of three si pseudofermion holes: 

(i') and (ii') transitions with the same values for the deviations AA"co, AMc^_i/2, 
AM^ _i/2, AMc, AMg as for the above general (i) and (ii) transitions, respectively, 
and AA^j\ = 3. The minimal excitation energy for such transitions is zero and 2/x, 
respectively. 
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For the Mott-Hubbard insulator quantum phase such that n = 1 one often shifts 
the ground-state zero-energy level to the middle of the Mott-Hubbard gap. Then the 
excitation energies 0, Emh, and 2Emh become 0, Emh/'^, and 3Emh/'2', respectively. 
Here we call three-holon states the excited energy eigenstates of the transitions (iii) 
because these involve the creation of three holons. Such excited energy eigenstates also 
involve the annihilation of three cO pseudofermions. Moreover, we call three-sl-hole 
states the excited energy eigenstates of transitions (i') and (ii'). These involve creation 
of three si pseudofermion holes. Finally, LHB and UHB excited energy eigenstates 
belonging to the sub-class of the general states generated by both the transitions (i) 
and (ii) such that AA^^^^ = 1 are called one-holon and one-sl-hole states. 

We start by evaluating the relative weight of the excited states of type (iii) (or states 
of types (iii), (i'), and (ii')) relative to the weight of the states of types (i) and (ii) (or 
states of types (i) and (ii) such that AA^^'\ = 1) when all these states are generated by 
application onto the ground state of the spin-up electron creation operator. (The types 
of excited energy eigenstates refer to the above transitions which generate these states 
from the ground state.) To assess the importance for finite values of L of the three-holon 
final states of type (iii) and three-sl-hole states of types (i') and (ii'), we perform an 
exact diagonalization of small chains. There is no one-to-one correspondence between 
the small-chain weight associated with each specific excited energy eigenstate and the 
excited-energy-eigenstate weight for L —>■ oo. However, there is such a correspondence 
for the spectral-weight sum rules of the states of types (i) and (ii) generated by dominant 
processes, relative to the sum rule of the states of type (iii). Moreover, within the 
excited energy eigenstates generated by dominant processes, we consider the relative 
weight of the states of types (i') and (ii'). For these spectral-weight sum rules the small- 
chain results provide values for the relative weights which agree up to 99% with the 
corresponding thermodynamic-limit values. 

The full electron addition and removal spectrum for six sites with six electrons (half 
filling) is shown in Fig. ^for U = 12t. (We checked that similar results are obtained for 
larger finite systems.) The first Hubbard bands at -^Emh/'^ are generated by dominant 
processes, whereas the second Hubbard bands at +3Emh/'^ and —3Emh/'2 result form 
processes generated by the above i > operators Ojj^^ and Oij^a, respectively. Note that 
the first and second Hubbard bands are well separated. The weights of the latter bands 
are orders of magnitude smaller than the contribution from the first Hubbard bands. 
The states centered around 3Emh/'^ are three-holon states of type (iii). As a result of 
the half-filling particle-hole symmetry, there is a corresponding structure for electron 
removal centered around —3Emh/'^- The latter structure is associated with creation of 
two rotated-electron unoccupied sites. (We recall that the —1/2 holon number deviation 
selection rule Q refers to excited-energy-eigenstate electronic densities such that n <1. 
There is a corresponding +1/2 holon number deviation selection rule for n> 1.) 

In figure |21 we plot the contribution of different excited energy eigenstates to the 
sum rule for half filling. For that, we have followed adiabatically the weights of different 
states as the value of U/ 1 is reduced, and summed the weight over the particular family 
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Figure 1. Electron addition {u > 0) and removal {uj < 0) spectrum for a half- filled 
finite-site ring, with U — 12 1. Note the logarithmic scale for the weights. The five- 
holon states are not shown. Indeed, the contributions of these states are extremely 
small and their energies are out of the shown energy window. 

of states. Analysis of the figure reveals that the contribution of the three-holon excited 
energy eigenstates of type (iii) to the total sum rule is largest at intermediate values 
U ~ 4t. It does not exceed 0.25% in the total sum rule. For large values of U/t it 
decreases as {t/U)'^, while for small values as {U/t)'^. The three-sl-hole contribution 
of the above states of type (ii') is also less than 0.25 %, and for small values of U/t it 
decreases as (U/t)'^. (These states belong to a sub-class of the excited energy eigenstates 
of type (ii).) While figure El corresponds to = 1, the relative spectral weight of the 
excited energy eigenstates of type (iii) decreases for decreasing values of the electronic 
density n. For instance, at quarter filling such a weight is 2% of that of half filling and 
vanishes as n — 0. Thus, at quarter filling and U ^ At the excited energy eigenstates 
of types (i) and (ii) (including the states of types (i') and (ii')) correspond to ~ 99.99% 
of the total spectral weight and the states (iii) to ~ 0.005% of such a weight. 

Evaluation for the model metallic phase of the available spectral- weight 
contributions by the method of Refs. [HI IH], confirms that the above results remain 
valid for L ^ oo. The exception is the relative weight of the excited energy eigenstates 
of types (i') and (ii'), which increases for increasing values of L. This is confirmed 
by the values given in the Table, which displays the relative weights generated by the 
transitions to the one-sl-pseudfermion-hole excited energy eigenstates that obey the 
relation (jHl) both for the one-electron removal and addition spectral functions. These 
weights were obtained for the one-electron spectral functions as U/t ^ oo. For L — >■ cxd 
and U/t ^ oo the use of the method of Ref . [12] leads to values of the relative weights 
for one-electron removal and addition of approximately 98% and 94%, respectively. The 
numbers provided in the Table confirm that also for smaller values of L the deviation 
value restrictions of Eq. ^ refer to a sub-class of excited energy eigenstates associated 
a substantial part of the full one-rotated-electron spectral weight. On the other hand. 



U = 12t 
n=1 
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Figure 2. The contribution of different states to the electron-addition sum rule for half 
filling. Over 99% of the sum rule is exhausted by one-holon and one-sl pseudofermion 
hole excited energy eigenstates of type (ii) . For larger systems this remains true if we 
also consider the states of type (ii'), which as the states of type (ii), are also associated 
with dominant processes generated by the i = operator ^. 



the set of excited energy eigenstates which obey the exact selection rule (H)) is larger: 
it corresponds to the whole one-rotated-electron spectral weight. The excited energy 
eigenstates of types (i') and (ii') are associated with dominant processes generated by 
the i = operator Oo,j,cr- contrast, the states of type (iii) are not generated by 
dominant processes and remain having very little spectral weight as L — oo. 

Let us next consider the limit of large U/t, where the L oo problem can 
be handled analytically by the method of Ref. ^21- In that case the general cO 
pseudofermions of Refs. [HI E] become the spin-less fermions used in the studies of 
Refs. [m I12j. Our first goal is the confirmation that the contribution from the three- 
holon states of type (iii) is very small and leads nearly to the same relative weight both 
for L ^ oo and finite values of L. For large values of U/t one can derive a systematic 
t/U expansion for the electron - rotated-electron unitary operator defined by Eqs. (21)- 
(23) of Ref. P^. (See reference which studies that transformation for large values 
of U/t.) By use of the inverse of the relation between the one-electron and one-rotated- 
electron operators given in Eq. (19) of Ref. jTHI in such a t/U expansion, we find after 
some algebra the following partial sum rule, 

J ^2UHB(^) duj=^(^^- 2SoSi - 2S1S2 - 2S0S2) 

X {nxo,conxi,conx2,co) , (7) 
where fix- co = fl cq fx co- The expectation value (| — 2SoSi — 2S1S2 — 2S0S2) 
refers to the spin degrees of freedom. The partial sum rule ((Tj) corresponds to the 
second UHB. This band is generated by transitions to the above three-holon states of 
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type (iii) and involves annihilation of three cO pseudofermions. The expectation value 
to find three neighboring local cO pseudofermions is, 

{fxo,cofxo,cQ} {fxo,cQ fxi,co} {fxo,cofx2,co} 
{fxi,cofxo,co) {fxi,cofxi,co) {fxi,cofx2,co) 
{ fx2 , cO fxo , cO ) ( fx2 , cO fxi , cO ) ( fx2 , cO fx2 , cO ) 



{^XQ, cO '^xi, cO ^X2, co) 



3 2nsin^(7m) ^ sin^(7rn) sin(27rn) 



^2 

n sin^(27m) 

This expectation value has the following limiting behavior, 



(8) 



if n <^ 1 : 



(^xo.cO ^xi,cO ^X2,co) — { (9) 

1 - 3 (1 - n) if 1 - n < 1 . 

Note that the spectral weight generated by the transitions to the three-holon states of 
type (iii) decreases rapidly away from half filling. At quarter filling (n = 1/2) it is about 
2% of that at half filhng. 

The expectation values associated with the spin degrees of freedom can in the 
thermodynamic limit and for large values of U/t be evaluated by use of the relation 
between the ID Hubbard model and the spin 1/2 isotropic Heisenberg chain |1H I12j. 
This leads to the following values 



so that, 



(5051) spin = ^ - In 2 ^ -0.443147 (10) 

(5052) spin = ^ - 4 In 2 + ^ C(3) ^ 0.182039 , (11) 



^-2SoSi- 28182-28082) =121n2-^C(3)^2.91. (12) 

^ / spin 



The sum rule of the second UHB is then given by, 



J ^™(cu)c?cu ^ 2.91 — {n,,, CO fix,, CO nx2,co) ■ (13) 

For the six-site finite-size cluster, the expectation value ()12|) is, 
3 



28081-28182-28082) = (169 + 17V13)/78 ^ 2.95, (14) 

^ I spin 

which is about 1% off from the thermodynamic-limit value given in Eq. p2|) . The 
asymptotic behavior 2.95 t^/t/'' is shown in Fig. |2las a dashed line. 

Our above numerical results for all values of V j t and a small system lead to the same 
general picture as the results for L 00. For electron addition the relative spectral 
weight of the excited energy eigenstates of types (i) and (ii) generated by dominant 
processes is minimum for U ~ 4t. This minimum value decreases with decreasing 
density. For half filling it is given by ~ 99.75%, whereas for quarter filling it reads 
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N rotated-electron removal rotated-electron addition 



6 


0.998792 


0.977515 


8 


0.997486 


0.972141 


10 


0.996277 


0.968088 


12 


0.995178 


0.964847 


14 


0.994176 


0.962156 


16 


0.993258 


0.959862 


18 


0.992409 


0.957867 


20 


0.991622 


0.956105 


22 


0.990886 


0.954531 


24 


0.990196 


0.953109 



Table 1. The relative weight of the one si pseudfermion hole contributions that obey 
relation © in the U/t oo limit for finite-size systems of N electrons. 



~ 99.99% and in the limit of vanishing density it becomes ~ 100.00%. The extremely 
small amount of missing one-electron addition spectral weight corresponds mainly to the 
three-holon states of type (iii) generated by the i > operators djj^a- Transitions from 
the ground state to higher-order five-holon/five-c-hole states generated by the latter 
operators lead to nearly vanishing spectral weight. 

The three-sl-hole states of types (i') and (ii') are generated by a sub-class of 
dominant processes. The relative weight of these excited energy eigenstates increases 
for increasing values of the system length L. For one-electron addition its maximum 
value occurs at half filling and is of about 6% as L ^ oo. (At n = 1 there are no LHB 
one-electron addition excited energy eigenstates of types (i) and (i').) For half filling, 
all values of U/t, and L ^ oo over 99% of the one-electron addition spectral weight 
corresponds to generation of one holon and one and three si pseudofermion holes. This 
is similar to the relative weights of Fig. |21 for L finite. For L —>■ oo the amount of 
one-electron spectral weight generated by dominant processes increases for decreasing 
values of the electronic density for all values of U/t. 

While for electron addition the states of types (i) and (ii) such that AA^^^^ = 1 
correspond to at least 94% of the spectral weight, for electron removal the AA^^^^ = 1 
excited energy eigenstates amount for at least 98% of the total weight. However, we note 
that concerning the small amount of spectral weight generated by the excited energy 
eigenstates such that AA^^'\ > 1, the case of the ID Hubbard model is different from 
that of the t — JxY model considered in Ref. [21]. The significant difference is the SU{2) 
symmetry in the spin sector of the ID Hubbard model, which by standard selection rules 
prohibits matrix elements that are present in the t — Jxy case. For instance, for low 
energy such a symmetry protects the lower-edge of the removal and addition spectrum of 
being dressed by si pseudofermion particle- hole excited energy eigenstates in the main 
conformal tower (see Fig. 5 in Ref. The low-energy AArj\ = 3 excited energy 

eigenstates are for the ID Hubbard model mostly associated with the next conformal 
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tower centered at 3k p for the spin excitations. In the t — Jxy case, there is no similar 
protecting mechanism and, therefore, the weight is more easily redistributed among the 
si pseudofermion particle-hole excited energy eigenstates. Let us also note that the 
studies of Ref. |22i have addressed the same question for other models. The authors 
of that reference found that the weight coming from AN^^ = 1 excited states is also 
dominant. 

3. ROTATED-ELECTRON GENERATORS OF THE ONE-ELECTRON 
SPECTRAL WEIGHT 

In this section we study the pseudofermion microscopic processes associated with the 
excitations cj^^IGS*) and Cj^(j\GS) which could be expressed in terms of operators of 
general form given in Eq. Q. Fortunately, for all values of U/t more than 99% of the 
spectral weight associated with such excitations corresponds to the i = contributions 
of expression Q. Therefore, here we limit our study to the i = operators Olj^a^ 
and Oo j^cr , which generate nearly the whole one-electron spectral weight measured 
in photoemission experiments ^ Indeed, for the electronic densities of the TTF 
and TCNQ stacks of molecules considered in the photoemission experiments of Ref. 
[Y] the studies of the previous section reveal that the dominant processes generated 
by the operators fi'J and 0o,j,a account for more than ^ 99.9% of the total one- 
electron spectral weight. We note that in spite of the recent improvements in the 
resolution of photoemission experiments [2], it is difficult to measure the finest details 
of the electronic structure experimentally, in part due to the extrinsic losses that occur 
on very anisotropic conducting solids such as the organic compound TTF-TCNQ |23j . 
Therefore, the less of 0.01% of missed theoretical spectral weight is irrelevant for the 
description of the spectral features measured by photoemission experiments. 

Following the exact selection rule of Eq. (0) and the corresponding rule for L^, _i/2, 
the CPHS ensemble subspaces of the excited energy eigenstates generated by application 
of the operators c]^ and Cj^o- onto the ground state can have the following values for 
the numbers {Lc, _i/2, Ls,-i/2}, 

c] I operator: {0, 0} , {1, 0} , {0, 1} , {1, 1} ; Cjj operator: {0, 0} ; 

I operator: {0, 0} , {1, 0} ; Cjj operator: {0, 0} , {0, 1} . 

We emphasize that the expressions of the operators Oqj^^ and Olj ^r o^^ly depend 
on the values {L^ -1/2, -^^s, -1/2} of the CPHS ensemble subspace they refer to. Different 
CPHS ensemble subspaces with the same values for {Lc, _i/2, -1/2} have the same 
expressions for the operators ^o,i,o- and 6o,j,cr- Evaluation of the commutators given in 
Eq. (27) of Ref. ^ for the one-electron case considered here together with the property 
that the i = operators 60. a and j ^ have the same expressions in terms of rotated- 
electron creation and annihilation operators as the corresponding operators Oj^cr and 
dj ^, respectively, in terms of creation and annihilation electronic operators, leads to the 
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following expressions, 



and 



<.,i-7==f^9-T, {1,0}, 

el I = ^ 5] . , {0, 1} , 

01-,= -^ ^ ' = g, I , {1,1}, (15) 



^o,jA = ^hU {0,0}, 



kl.l=i,.l, {0.0}, (16) 

where we recall that < < Na and < < N^. Here the values of the 
numbers {L^ _i/2, -^^5,-1/2} are those provided above and A^°, N^, and N'^ are the 
electron numbers of the initial ground state. 

The CPHS ensemble subspaces associated with ?7-spin value deviations ASc and spin 
value deviations ASg such that ASc = and ASs = —1/2 are not permitted for 

initial ground states such that A^*^ = Na and = N^, respectively. The reason is that 
for such ground states the ?7-spin and spin values are 5"° = and S^ = 0, respectively, 
and thus negative deviations ASc and ASg are not allowed. It follows that for = 
ground states, the transitions to CPHS ensemble subspaces such that Ls,_i/2 = cannot 
be generated by application of the operators ^ and Cj_ | onto such initial ground states. 
The case of the iV| = N'^ ground states plays an important role in the applications of 
our results and those of Refs. |Hl IH]- In the following, we provide the expression of 
the operators ^0,^,0- and Oqj^^ of Eqs. (fT5|l and (fT6|) in terms of pseudofermion creation 
and annihilation operators for the set of CPHS ensemble subspaces associated with the 
= ground states. 

The one-electron processes generated by the rotated-electron operator c] give rise 
to the LHB when Mc^ -1/2 = and UHB when Mc^-1/2 = 1- According to the exact 
selection rules (jH), those are the only permitted values for the excited-energy-eigenstate 
—1/2 holon numbers. Since the numbers of the initial ground-state CPHS subspace 
are known and well defined, here we characterize the CPHS ensemble subspaces of the 
excited states associated with the dominant processes by the deviation numbers and 
numbers AA^^co, AA^^i, ANs2, and //s, -1/2- For the UHB we also consider the numbers 
AA'ci and Lc, -1/2- As discussed above, our aim is the study of the dominant microscopic 
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processes that generate the excitations Cj ^\GS) and Cj^i\GS) where \GS) is a = 
ground state. 

For simphcity, here and in the ensuing section we consider local generators which 
are the Fourier transform of the corresponding generators associated with processes such 
that for each CPHS ensemble subspace all pseudofermion and pseudofermion holes are 
created away from the Fermi points for the cO and si branches and from the momentum 
values of largest absolute value for the s2 and cl branches. Similar expressions can be 
derived for processes including creation of pseudofermions or holes at such Fermi points 
and limiting momentum values [Hlini- 

The whole Mc,_i/2 = LHB spectral weight of the excitation ICS") corresponds 
to {Lc, -1/2, -^5,-1/2} = {0, 0} CPHS ensemble subspaces. According to Eq. (fTBj) . in such 
subspaces the above excitation reads 9l j^\GS) = c] j (1 — nj^i)\GS) where (1 — n^, j.) 
is the LHB projector. We recall that the LHB processes do not exist for the = Na 
half-filling ground state. Otherwise, the spectral weight of this excitation is generated 
by transitions to the excited energy eigenstates which span the set of CPHS ensemble 
subspaces such that, 

AN,o = l, AN,i = -2N,2, Ns2 = 0,l. (17) 

For simplicity, here and in other equations given below we included only the au branches 
with finite pseudofermion occupancy for the subset of CPHS ensemble subspaces 
considered in our study, which correspond to nearly the whole spectral weight. Thus, 
the deviation AA^^'^ = AA^^co — 2AA^si — 2Ns2 in the number of si pseudofermion holes 
is AN^^ = — 1 + 2Ns2 for these subspaces. We note that for the present case of zero- 
magnetization the si pseudofermion band is full for the ground state. Therefore, the 
deviation AA^^^^ equals the number of si pseudofermion holes of the CPHS ensemble 
subspace excited energy eigenstates. 

At least 94% of the spectral weight of 9lj^-^\GS) = ^ {l—nj^i)\GS) corresponds to 
the CPHS ensemble subspace such that AA'^o = 1, ANsi = 0, Ns2 = 0, and AA'j\ = -1. 
These numbers obey the relation (jHI). In this subspace the pseudofermion expression of 
the operator O^ j^-^ = cj ^ (1 — fij |) is given in Eq. (47) of Ref. p. Nearly the whole of 
the remaining spectral weight of the above LHB excitation is associated with the CPHS 
ensemble subspace whose numbers are AA'co = !> ANgi = —2, Ns2 = 1, and AA'^'\ = 3. 
Such values do not obey the relation In this subspace, the pseudofermion expression 
of the operator O^ j ^ = ^ (1 — fij^i) reads, 

^l.,T = 4 T (1 - = e-^^^^^^ fl^,2 f.,,„si flco , (18) 

where the values of the phase-factor momentum APj and U/t independent real positive 
constant Gc are specific to the subspace and are given in Ref. the effective si lattice 
index reads j' = jn/2 and the effective s2 lattice is reduced to a single site such that 
Xi = Xj. Here and below the equality j' = j n/2 (and j' = j 6 for the effective cl lattice) 
should be understood as j' being the closest integer number to jn/2 (and j6). (We 
recall that the effective cO lattice index j equals that of the rotated-electron lattice |17j.) 
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Next, we consider the excitation Cj^i\GS). Again, more than 99% of the spectral 
weight corresponds to the excitation 6qj^^\GS). Such an excitation is associated with 
{-^c, -1/2, -^5,-1/2} = {0, 0} CPHS ensemble subspaces where according to Eq. (|T!)jl . 
9oj^l\GS) = Cj^i (1 — njj)\GS). Here the projector (1 — n^j) is associated with the 
ground-state lack of rotated-electron double occupation [T3j. Most of the spectral weight 
of this excitation is generated by transitions to the states which span the set of CPHS 
ensemble subspaces such that, 

AN,o = -l, ANsi = -1 - 2Ns2 , Ns2 = 0,l, (19) 

and AN^^ = l+2iV,2. Up to 98% of the spectral weight of9oj^i\GS) = c^- | (l-n^- 
corresponds to the CPHS ensemble subspace whose deviations read ANco = —I, 
ANsi = -1, A^^2 = 0, and AN^^ = 1. These values obey the relation ©• The 
pseudofermion expression of the operator Oqj^i = Cj,|(l — nj,-\) in this subspace, is 
given in Eq. (46) of Ref. jU]. Nearly the whole of the remaining spectral weight of 
the above electron-removal excitation is associated with the CPHS ensemble subspace 
whose numbers read ANco = —1, AA^^i = —3, Ns2 = 1, and AA^^^^ = 3. These do not 
obey the relation ®. In this subspace the pseudofermion expression of the operator 
^o,j,i = (1 - nj- t) is given by, 

where j' = j n/2 and Xi = Xj. 

For the Mc^-1/2 = 1 UHB excitations, let us start by considering electronic densities 
such that < 5 < 1, where 6 = [1 — n] is the doping concentration. We are assuming 
that 6 can be small, but not vanishing. The case of vanishing doping concentrations 
corresponds to the Mott-Hubbard insulator phase where {Na — N^) = and to the 
doped Mott-Hubbard insulator such that {N^ — N^) is finite and will be addressed 
in the ensuing section. For finite doping concentrations the whole UHB spectral 
weight of the excitation 6lj^^\GS) corresponds to {Lc, _i/2, i^s, -1/2} = {0, 0} CPHS 
ensemble subspaces. According to Eq. (fT6|l . such an excitation can be written as 
9lj^-^\GS) = Cj^-^hj^i\GS) where hj^i is the UHB projector. For finite values of S, most 
of the spectral weight of this excitation is generated by transitions to the excited energy 
eigenstates which span the set of CPHS ensemble subspaces such that, 

AAr,o = -l, ANsi = -1 - 2Ns2 , N,^ = 1 , Ns2 = 0,l- (21) 

The deviation in the number of si pseudofermion holes of these subspaces is AA^^'\ = 
1 + 2N,2- At least 94% of the spectral weight of the excitation Oq^j^-^IGS) — cj | 

\GS) 

corresponds to the CPHS ensemble subspace such that AA^co = —1, AA^^i = — 1, 
Aci = 1, Ns2 = 0, and AA"^^ = 1. These obey the relation (jH)). In this subspace, 
the pseudofermion expression of the operator O^j^^ = cj^ j is given by, 

^l.,T = 4t %i = e-^^'^''^ fl„ci /x„co 4„.i , (22) 
where j' = jn/2 and j" = j 6. Nearly the whole of the remaining spectral weight of the 
above UHB excitation is associated with the CPHS ensemble subspace whose numbers 
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are AiVco = —1, AA^^i = —3, = 1, = 1, and AN^^ = 3. These however do not 
obey the relation ®. In this subspace the pseudofermion expression of the operator 

^o,j,T = ^T^i.i is Si^^^ by, 

n 
2Gc 

where j' = j n/2, j" = j 6, and xi = Xj. 



4. THE METAL - MOTT-HUBBARD INSULATOR QUANTUM PHASE 
TRANSITION 

The further understanding of the microscopic processes behind the quantum phase 
transitions is a problem of great physical interest, as mentioned in Sec. I. For the ID 
Hubbard model there is only one quantum phase transition as a function of the on-site 
repulsion U, which occurs at f/ = for half filling. As a function of the interaction, such 
transitions are controlled by the dependence on that interaction of the local quantum 
entanglement ^3]. On the other hand, for all finite values of U there is a metal - 
insulator quantum phase transition which occurs as a function of the electronic density 
at n = 1. In terms of the doping that quantum phase transition occurs at 5 = 0. Here 
we study the microscopic effects of that quantum phase transition on the one-electron 
spectral properties. 

Let us consider the case when S is vanishing, i.e. there are no cO pseudofermion 
holes (half-filling) or there is a finite number of such holes in the initial ground 
state (doped Mott-Hubbard insulator). As the number of holes decreases and one 
reaches the doped Mot-Hubbard insulator regimen, there arises a competition of the 
UHB processes considered in the previous section with other UHB processes. The 
latter processes correspond to the excitation 9lj^.^\GS) for CPHS ensemble subspaces 
such that {Lc, -1/2, -Z^s, -1/2} = {1; 0}. Thus, according to Eq. (fTBj) . one has that 
elj^^lGS) = -[{-ly/VNa - iV" + 1] d]^^nj,^\GS), where n^- t is the UHB projector. 
Note that for finite values of the doping concentration 6 the one-electron spectral weight 
associated with such an excitation vanishes in the themodynamic limit L — > cxd. Most 
of that weight is generated by transitions to the excited energy eigenstates which span 
the set of CPHS ensemble subspaces such that, 

AiV,o = -l, AiV.i = -1 - 2iV,2 , L,,_i/2 = 1, iV,2 = 0,l, (24) 

and ANgi = 1 + 2Ns2- For vanishing values of the doping concentration at least 94% 
of the spectral weight of the excitation Olj^-^IGS) corresponds to the CPHS ensemble 
subspace with numbers ANco = —1, AiVsi = —1, _i/2 = 1, Ns2 = 0, and AN^^ = 1. 
Such numbers obey the relation ®. In this subspace the pseudofermion expression of 
the operator Ol reads, 

^o,.,T = -T^T^f^Su^T = -^^^_^,^^fx„sifx„co, (25) 



Spectral microscopic mechanisms and. 



17 



where j' = j/2. Nearly the whole of the remaining spectral weight of such a 
UHB excitation is associated with the CPHS ensemble subspace whose numbers are 
AA^co = —1) ^Nsi = —3, i^c, -1/2 = 1, and AA^j\ = 3. These do not obey the relation 
(jH)). In this subspace the pseudofermion expression of the operator 6^ ■ ^ is given by, 



- (-^)^ -^t ^.^ 



-ij[APj+n] ^ 

fxi,s2 fx-ij^2^sl fx-ij^^,sl fx-i,sl fxj,cO ; (26) 



2GcVNa -N^ + l 
where j' = j/2 and Xi = Xj. 

For finite values of 6, the UHB processes associated with CPHS ensemble 
subspaces such that {Lc^-1/2, -^^s, -1/2} = {Ij 0} are allowed. However, for finite 
doping concentrations the relative UHB weight of the {Lc, _i/2, -^^5,-1/2} = {0, 0} and 
{Lc, -1/2, Ls^ ^1/2} = {1, 0} processes is approximately given by 1/6 L and thus vanishes 
in the thermodynamic limit, as mentioned above. Therefore, the weight of the UHB 
{Lc, -1/2, -Ls, -1/2} = {1, 0} processes associated with the CPHS ensemble subspaces of 
numbers ()24|) vanishes for finite values of 6. 

In contrast, for half filling only the UHB {1, 0} processes contribute, whereas for 
finite yet small values of {Na — N^) both the {0, 0} and {1, 0} processes contribute to 
the UHB spectral weight. If the value of {Na — N^) further increases so that 6 becomes 
finite, the weight of the UHB {1, 0} processes vanishes. For 7^ Na and {Na — A^°) 
finite but small, the relative weight of the {0, 0} and {1, 0} processes is approximately 
given by 1/[A^„ - A^o + 1]. 

We thus conclude that the collapse of the LHB processes and the interplay between 
the {0, 0} and {1, 0} UHB processes for decreasing values of {Na — A^°) are the main 
effects of the (A^^ — N^) — >■ metal - Mott-Hubbard insulator quantum phase transition 
onto the one-electron spectral properties. The ground-state - excited-energy-eigenstate 
transitions associated with the {0, 0} and {1, 0} UHB processes change the value of rj 
spin by Ar/ = — 1/2 and Arj = +1/2, respectively. For the metallic phase at finite values 
of the doping concentration the whole UHB weight corresponds to excited states with 
deviation Arj = —1/2, whereas for the Mott-Hubbard insulator phase only the excited 
states with deviation Arj = +1/2 are allowed. However, the physics of the doped Mott- 
Hubbard insulator, which corresponds to a finite number of holes {Na — N'^), is different: 
It involves a competition between the two above classes of states. Such a competition 
is mainly controlled by the form of the operators given in Eqs. ()22j) and which 
generate excited energy eigenstates with ?7-spin deviations Arj = —1/2 and Arj = +1/2, 
respectively. 



5. CONCLUDING REMARKS 



In this paper we used pseudofermion description in the study of the non-perturbative 
microscopic processes that control the unusual one-electron spectral- weight distributions 
of the ID Hubbard model. Our results are useful for the further understanding and 
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description of the microscopic mechanisms behind the unusual finite-energy spectral 
properties observed in quasi-lD compounds ^1 |2j. They also provide new insights 
about the microscopic mechanisms of quantum phase transitions in electronic correlated 
problems. We found that for the doped Mott-Hubbard insulator there is a competition 
between the microscopic processes which determine the UHB one-electron spectral 
distributions of the Mott-Hubbard insulator phase and metallic phase for finite values 
of the doping concentration. 
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